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ABSTRACT 

We introduce and discuss the properties of a theoretical {B—K)-( J—K) integrated 
colour diagram for single-age, single-metallicity stellar populations. We show how this 
combination of integrated colours is able to largely disentangle the well known age- 
metallicity degeneracy when the age of the population is greater than ~ 300 Myr, 
and thus provides valuable estimates of both age and metallicity of unresolved stellar 
systems. We discuss in detail the effect on this colour-colour diagram of o-enhanced 
metal abundance ratios (typical of the oldest populations in the Galaxy), the presence 
of blue horizontal branch stars unaccounted for in the theoretical calibration, and of 
statistical colour fluctuations in low mass stellar systems. In the case of populations 
with multiple stellar generations, the luminosity-weighted mean age obtained from this 
diagram is shown to be heavily biased towards the youngest stellar components. We 
then apply this method to several datasets for which optical and near-IR photometry 
are available in the literature. We find that LMC and M31 clusters have colours which 
are consistent with the predictions of the models, but these do not provide a sensitive 
test due to the fluctuations which are predicted by our modelling of the Poisson 
statistics in such low-mass systems. For the two Local Group dwarf galaxies NGG185 
and NGG6822, the mean ages derived from the integrated colours are consistent with 
the star formation histories inferred independently from photometric observations of 
their resolved stellar populations. 

The methods developed here are applied to samples of nearby early-type galaxies 
with high quality aperture photometry in the literature. A sample of bright field and 
Virgo cluster elliptical galaxies is found to exhibit a range of luminosity-weighted 
mean ages from 3 to 14 Gyr, with a mean of ~8 Gyr, independent of environment, 
and mean metallicities at or just above the solar value. Colour gradients are found in 
all of the galaxies studied, in the sense that central regions are redder. Apart from 
two radio galaxies, where the extreme central colours are clearly driven by the AGN, 
and one galaxy which also shows a radial age gradient, these colour changes appear 
consistent with metallicity changes at a constant mean age. Finally, aperture data for 
five Virgo early-type dwarf galaxies shows that these galaxies appear to be shifted to 
lower mean metallicities and lower mean ages (range 1 to 6 Gyr) than their higher 
luminosity counterparts. 

Key words: Galaxies: clusters: individual (Virgo) - galaxies: stellar content - galaxies: 
elliptical and lenticular - galaxies: dwarf 


1 INTRODUCTION 

Within the currently favoured ACDM models of galaxy and 
large scale structure formation, dwarf galaxies play a crucial 
role. The first stars are predicted to form in small dark mat¬ 
ter halos of mass ~ 10® Mq at redshifts of a « 20 — 30 (see 


Bromm and Larson 2004). These small dark matter halos 
subsequently form the building blocks from which the larger 
galaxies we see around us today are assembled (Kauffman 
et al. 1993; Cole et al. 2000; Mathis et al. 2002). As with 
every building site many of the smaller component parts are 
left behind once the larger structures are formed and so the 
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ACDM model predicts that there should be many small dark 
matter halos around us today. 

If this debris of small dark matter halos can be associ¬ 
ated with dwarf galaxies then there is a large discrepancy 
between what is observed and what is predicted by the¬ 
ory, this is known as the sub-structure problem (Moore et 
al. 1999; Klypin et al. 1999). The standard ACDM model, 
without dwarf galaxy formation suppression processes, pre¬ 
dicts many more small dark matter halos (with a low mass 
slope of the mass function of a « — 2) compared to obser¬ 
vations of the faint end of the global luminosity function 
of galaxies (a « —1.2; Blanton et al. 2001, Norberg et al., 
2002). There have been numerous suggestions for a solution 
to this problem. These typically involve the prevention of 
gas falling into small dark matter halos (Efstathiou 1992), 
the removal of gas via supernova winds from the hrst gen¬ 
eration of stars (Dekel and Silk 1986) or the ‘squelching’ 
of dwarf galaxy formation by ionising photons (Tully et al. 
2002). These mechanisms prevent or severely inhibit star for¬ 
mation in small halos; the halos still exist, but they remain 
hidden from us. Interest in the possiblity that star formation 
in dwarf galaxies has in some way been delayed compared 
to larger galaxies has come from the extensive studies made 
using Sloan Digital Sky Survey data. This seems to indi¬ 
cate that it is the low luminosity galaxies that contain the 
youngest stars (Kauffmann et ah, 2003). It is difficult to see 
how this ‘downsizing’ of the galaxy population hts in with 
the standard ACDM model (Kodama et al. 2004), but it is 
an observation we can test with the methods we discuss in 
this paper. Recently other solutions have been suggested, 
for example by adjusting the initial dark matter fluctuation 
spectrum, something suggested by combining the WMAP, 
2dF and Lyman-a forest observations (Spergel et al. 2003), 
small dark matter halos themselves become rare. 

By studying in detail the sub-structure problem in dif¬ 
ferent galaxy environments we can hope to gain an insight 
into the mechanisms that prevent dark matter halos reveal¬ 
ing themselves to us as dwarf galaxies or comment on how 
universally rare small halos might be. We have previously 
quantffied the relative numbers of dwarf galaxies in differ¬ 
ent environments and find many more dwarf galaxies in the 
Virgo and Fornax clusters than can be found in the gen¬ 
eral field (Phillipps et al. 1998a; Sabatini et al. 2003, 2004; 
Roberts et al. 2004, 2005; Davies et al. 1988; Kambas et al. 
2000). The general field result is also consistent with obser¬ 
vations of the dwarf galaxy population of the Local Group 
(Mateo 1998; Pritchet & van den Bergh 1999). It appears 
that global star formation suppression mechanisms do not 
work. The environment plays a major role in either pro¬ 
ducing a large number of additional dwarf galaxies or in 
revealing more of the existing dark matter halos to us. 

The expected properties of the first galaxies to form 
within the ACDM model seem to be more closely matched 
to the properties of globular clusters (a star formation event 
followed by almost total gas expulsion) than it does to the 
more complicated star formation histories of dwarf galax¬ 
ies. Although dwarf galaxies appear to be relatively simple 
stellar systems (there is little evidence for spiral stellar den¬ 
sity waves, bars or jets, but see Barazza et al. 2002) they 
can be more susceptible to exterior influences. In the main 
dwarf galaxies show a morphology density relation with the 
gas poor dSph galaxies being closer to the bright galaxies in 


the Local Group and closer to the centres of galaxy clusters 
than are the dirr galaxies (Ferguson & Sandage 1989; Mateo 
1998; Boyce et al. 2001; Sabatini et al. 2003; Roberts et al. 
2004). Within a cluster environment, or when they are close 
to a more massive galaxy, dwarf galaxies may be subject 
to ram pressure and/or tidal stripping, tidal destruction or 
‘harassment’ processes (for a review see Sabatini et al. 2004). 

Currently we essentially only have detailed star forma¬ 
tion histories for Local Group dwarf galaxies because these 
are the only ones for which we have resolved observations 
of their stellar populations. These have been interpreted us¬ 
ing stellar colour magnitude diagrams and models that im¬ 
plement different stellar ages and metallicities (Grebel & 
Stetson 1999; Smecker-Hane et al. 1996; Gallart et al. 1999; 
Mighell and Rich 1996; Skillman et al. 2003). Each of the 
Local Group dwarf galaxies has a quite different star forma¬ 
tion history when studied in detail, with stellar populations 
covering a wide range of ages and metallicities. In general 
the dSph galaxies have no young stars; the stellar population 
is dominated by stars that are older than 10 Gyr, but there 
are some intermediate age stars (TIO Gyr) indicating that 
star formation continued for many Gyr after their forma¬ 
tion (Grebel 2002). The current, most widely used, AGDM 
solution to the substructure problem (that small galaxies in 
the early Universe should experience a single catastrophic 
episode of star formation) is not supported by observations 
of the star formation histories of Local Group dwarf galaxies 
or by the large numbers of dwarf galaxies found in clusters. 

Although subject to similar physical processes dwarf 
galaxies in clusters may have very different origins and his¬ 
tories to those in the Local Group. Within the cluster en¬ 
vironment there are many more dwarf galaxies per giant 
(Phillipps et al. 1998b, Sabatini et al. 2003), so the clusters 
are not being assembled from groups like the Local Group. 
Also unlike the Local Group there are many early type (dE) 
dwarf galaxies that are not companions of the giant galax¬ 
ies (Sabatini et al. 2003, 2004; Roberts et al. 2004). Ide¬ 
ally we would like to study the star formation histories of 
these galaxies in the same detail as has been done for the 
Local Group galaxies to see if they are consistent with each 
other and with the properties predicted by AGDM - predom¬ 
inantly old, gas poor, low metallicity galaxies. Currently this 
is not possible for dwarf galaxies in the two nearest clusters 
to us (Virgo and Fornax) for two reasons. Firstly, because 
the bulk of the stellar population cannot be resolved at these 
distances (although the resolution of galaxies at these dis¬ 
tances is a key science driver for the proposed Extremely 
Large Telescopes). Secondly, many of the dwarf galaxies we 
are particularly interested in are of such low surface bright¬ 
ness (LSB) that it is impossible to obtain spectra and de¬ 
rive line indices (Worthey 1994, Vazdekis and Arimoto 1999, 
Kuntschner 2000). Thus, at present we are forced to use 
population indicators that can be applied to completely un¬ 
resolved LSB stellar populations. 

In the present study we explore the use of optical/near- 
IR colours as diagnostics of luminosity-weighted mean ages 
and metallicities of stellar populations. We demonstrate that 
such colours have similar diagnostic power to spectral line 
indices, while being much less demanding in terms of ob¬ 
serving time, data reduction and analysis. The development 
of such indicators is particularly timely with the recent re¬ 
lease of the full Two-Micron All-Sky Survey (Jarrett et al. 
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2000; henceforth 2MASS), and the prospect of deeper sur¬ 
veys such as the United Kingdom Infrared Deep Sky Survey 
(UKIDSS) recently commenced with the United Kingdom 
Infrared Telescope (UKIRT). 

Model predictions of optical/near-IR colours have previ¬ 
ously been presented by several groups (e.g. Fioc and Rocca- 
Volmerange 1997; Small et al. 2001, Girardi et al. 2002; 
Bruzual and Chariot 2003) and observational studies have 
already highlighted the potential of such measurements for 
breaking the age-metallicity degeneracy, which is the first re¬ 
quirement of any successful population diagnostic method. 
Peletier, Valentijn & Jameson (1990) used (U — V), (B — V) 
and (V — K) colours to constrain the ages and metallici- 
ties of 12 elliptical galaxies, while Peletier & Balcells (1996) 
demonstrated that (U — R) — (R — K) colour-colour plots 
gave useful information on the mean ages of galaxy bulges 
and disks. Bell et al. (2000) and Bell & de Jong (2000) 
used {B — R) and (R — K) colours in the analysis of the 
star formation histories of low surface brightness and spi¬ 
ral galaxies respectively. Puzia et al. (2002) and Hempel & 
Kissler-Patig (2004) have recently demonstrated the power 
of iy — I) — {y — K) colour-colour diagrams applied to analy¬ 
sis of the star-formation histories of unresolved extragalactic 
globular clusters. 

In this paper we first develop our own predictions of 
the optical/near-IR colours of stellar populations as a func¬ 
tion of age and metallicity, using models developed by two 
of the authors (MS and SC). We identify the {B — K)- 
(J — K) colour-colour plane as a particularly sensitive one 
for breaking age-metallicity degeneracy, and use photome¬ 
try to test these predictions for systems (globular clusters 
and two dwarf galaxies in the Local Group) with existing 
independent measurements of age and metallicity. We also 
apply these methods to a small number of Virgo cluster el¬ 
liptical and dwarf elliptical galaxies which have sufficiently 
high quality photometry in the 2MASS catalogues and B- 
band photometry in matched apertures, in the first stage in 
a study of the star formation history of cluster galaxies. 


2 THE THEORETICAL CALIBRATION 

We employed the stellar evolution model library by Pietrin- 
ferni et al. (2004), that spans the metallicity (scaled solar 
metal distribution) range between Z=0.0001 and Z=0.04 
(10 different metallicities in total); this model grid enables 
one to compute isochrones for ages from ~ 20-40 Myr up 
to 15-20 Gyr, covering all evolutionary phases from the 
Zero Age Main Sequence (ZAMS) up to the beginning of 
the Thermal Pulse (TP) phase along the Asymptotic Giant 
Branch (AGB) or carbon ignition, depending on the mass 
of the most evolved objects. In order to compute correctly 
the near-IR integrated colours of stellar populations older 
than ~100 Myr we have extended the evolutionary tracks 
of the masses that evolve along the AGB phase, until the 
end of the TP phase. This has been performed by means 
of the synthetic AGB evolution technique (see e.g. Marigo 
et al. 1996 and references therein) starting from the stel¬ 
lar structures at the beginning of the TPs. The use of this 
technique suffices for our purposes and avoids extremely 
time-consuming evolutionary computations. The analytical 
relationships necessary for the synthetic AGB computations 


(e.g. an equation for the rate of growth with time of the 
carbon-oxygen core, the core mass-luminosity relationship, 
etc) are taken from Wagenhuber & Groenewegen (1998); 
they are obtained from full evolutionary AGB calculations 
covering the relevant range of initial masses and metallic¬ 
ities. The analytical formulae for the determination of the 
model Teff along the extended AGB phase are not given ex¬ 
plicitly in Wagenhuber & Groenewegen (1998), but can be 
found in Wagenhuber (1996); the mass loss rate along the 
TP phase is computed according to the methods of Vas- 
siliadis & Wood (1993). The synthetic AGB evolution is 
stopped when the models start to cross the Hertzsprung- 
Russell diagram towards their white dwarf cooling sequence. 

Using this extended set of evolutionary tracks we com¬ 
puted isochrones for ages running from 40 Myr to 14 Gyr, 
suitably transformed to various photometric filters using the 
bolometric corrections and colour transformations described 
in Pietrinferni et al. (2004). For the AGB TP phase we have 
used the transformations described in Girardi et al. (2002) 
modified in the carbon star regime to reproduce the empiri¬ 
cal {B — K) — {J — K) colour-colour relationships for Galactic 
and LMC clusters’ carbon stars, as obtained from data by 
Mendoza & Johnson (1965) and Frogel & Cohen (1982). 

Integrated magnitudes and colours for the simple stellar 
populations (SSPs) described by our isochrones have been 
determined using a Salpeter (1955) Initial Mass Function 
(IMF). We adopted as lower and upper IMF mass cutoffs 
the values O.IMq and lOOM©, respectively, as in Bruzual 
& Chariot (2003). The choice of these cutoffs is important 
for the determination of the total mass of the system and 
the integrated magnitudes, but it is not critical for the inte¬ 
grated colour-colour calibration discussed in this paper, due 
to the following two reasons. First, the integrated colours 
are independent of the total stellar mass, apart from the 
case of stellar populations with masses typical of star clus¬ 
ters (see next subsection). Only in this case does the value 
of the total mass of the stellar population play a role in the 
determination of the integrated colours. Second, the upper 
limit of the mass contributing to the total flux and colours 
is fixed by the most massive star evolving at the age of the 
population, and it is at most around IOMq in our discus¬ 
sion. As for the lower mass limit, we have found that in the 
photometric filters of interest, MS objects less massive than 
~O.5M0 do not provide an appreciable contribution to the 
integrated flux. To give a quantitative estimate of the con¬ 
tribution of these low-mass stars, we have recomputed the 
(B — K) and (J — K) colours neglecting stars with masses be¬ 
low 0.5Mq . We found that the (J — K) values are decreased 
by less than 0.01 mag and {B — K) is decreased by only 
0.03-0.05 mag over the full range of ages and metallicities of 
our calibration. 

Figure 1 displays the theoretical integrated (B — K) and 
(J — K) colours for a subset of our isochrones with ages of, 
respectively, 40, 300 and 600 Myr, and 1, 3, 6, 10, and 14 Gyr 
(from bottom to top) and (from left to right) [Fe/H]=-1.79, 
-1.27, -0.66, 0.06, and 0.40. This combination of colours 
is able - for ages larger than ~ 600 Myr - to appreciably 
disentangle the effect of age and metallicity on the observed 
integrated fluxes. Some residual degeneracy is still present 
between 600 Myr and 1 Gyr at the high metallicity end of 
the calibration. The integrated colour (J — K) appears to 
be a good tracer of the initial metallicity of a Simple Stellar 
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Figure 1. Theoretical integrated {B — K) and (J — K) colours 
for isochrones with ages of 40 Myr, 300, and 600 Myr,l, 3, 6, 10 
and 14 Gyr (open circles), and metallicities [Fe/H] from —1.79 
(leftmost solid line) to +0.40 (rightmost solid line). Dashed lines 
connect points corresponding to SSPs of constant age. The arrow 
indicates the direction of the reddening vector. 


Population (SSP) and is weakly affected by age, whereas the 
integrated {B — K) is a good age indicator, mildly affected 
by the SSP metallicity. The direction of the reddening vector 
is also shown in the figure. 

We have studied the contribution of the individual evo¬ 
lutionary phases to the B, J and K integrated magnitudes 
and found, not unexpectedly, that the J and K integrated 
fluxes are dominated by AGB stars when the SSP age is 
below ~ 1 Gyr, and by upper Red Giant Branch (RGB) ob¬ 
jects for higher ages. As for the integrated B flux, the main 
contribution comes always from the upper Main Sequence 
(MS) and TO stars. This means that the integrated {J — K) 
colour is mainly determined by the colour of AGB and/or 
RGB stars, whose location is strongly affected by their initial 
metallicity, whereas {B — K) is sensitive to the magnitude 
and colour of the TO, hence to the SSP age. 

It is also very interesting to notice how the theoretical 
{B — K) — {J — K) integrated colour diagram closely mirrors 
the behaviour of the theoretical calibrations of pairs of Lick 
indices (e.g. Worthey 1994) like, for example Hjs — Fe5270, 
where Hp is a tracer of the stellar population age and 
T'e5270 is sensitive mainly to the metal content. 

For ages above ~10 Gyr the {B — K) colour tends to 
lose sensitivity to age, whereas below ^300 Myr there is full 
degeneracy between age and metallicity, since the lines cor¬ 
responding to different metallicities overlap. These general 
properties of the integrated [B — K)-{J — K) diagram are 
not peculiar to our set of models; they are displayed also i.e., 
by the Girardi et al. (2002) and Bruzual & Gharlot (2003) 
integrated colours. As an example, we show in Fig. 2 the 
Girardi et al. (2002) theoretical diagram, for the age range 
where the age-metallicity degeneracy is broken. 



Figure 2. Theoretical integrated (B — K) and (J — K) colours 
from Girardi et al. (2002) models. The calibration is shown for 
the age range where the age-metallicity degeneracy is broken. 

2.1 The effect of statistical fluctuations 

Given that the fast evolving upper RGB and AGB stars are 
the main contributors to the (J — K) integrated colours, 
it is very important to assess the effect of small number 
statistics on the integrated {B — K) and (J — K) colours. 
When the mass of the SSP under scrutiny is small - total 
masses up to 10® — 10® M©, the upper end of the mass 
spectrum of star clusters - the number of stars in these fast 
evolutionary phases is subject to sizable fluctuations from 
one SSP to another, with potentially large fluctuations of the 
integrated magnitudes and colours. To quantify this effect 
we have performed a series of Monte-Garlo simulations. 

We have considered the following combinations of ages 
(t) and [Fe/H] values: [Fe/H]=—1.79 and t =10 Gyr, 
[Fe/H]=—0.66 and t =10 Gyr, [Fe/H]=—1.27 and 
t =600 Myr. For each one of these combinations, we have 
drawn stars randomly from a Salpeter IMF and placed them 
in their evolutionary phases along the isochrone at that age, 
until a cluster mass of 10® M© is reached. The integrated 
colours of the resulting populations have been then com¬ 
puted. In general, we found that the fluctuations are negli¬ 
gible in B since this wavelength range is dominated by the 
much more populous MS phase. 

The open squares, triangles and pentagons in Fig. 3 
display the integrated colours obtained from 100 realizations 
each for the three adopted age and metallicity combinations. 
The (J — K) colours show a very large spread due to statisti¬ 
cal fluctuations of the number of stars along the upper RGB 
and AGB phases. This spread increases for decreasing age 
because of the shorter timescales. Hence, there are larger 
number fluctuations along the AGB phase that dominate 
the (J — K) colours at 600 Myr. The colour spread causes 
a 3 (t uncertainty of «2 dex in the inferred metallicity at 
600 Myr, and « Idex at 10 Gyr. The fluctuation of {B — K) 
is entirely due to the fluctuation of the K magnitudes and 
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Figure 3. Model colours with overlaid points showing the results 
of 100 simulations of 10® Mq synthetic stellar populations. The 
scatter in the colours of the synthetic populations stems from 
the statistical fluctuations in the numbers of the most luminous 
stars. Open squares, triangles and pentagons display populations 
with, respectively: [Fe/H]=—1.79 and t=10 Gyr, [Fe/H]=—0.66 
and t=10 Gyr, [Fe/H]=—1.27 and t=600 Myr. 



Figure 4. Observed colours of LMC clusters (points) overlaid 
on our model predictions for [Fe/H] between —1.79 and +0.06. 
Different point types denotes the mean cluster ages as determined 
by van den Bergh (1981). The effect of an additional extinction 
Ay=0.5 mag is also displayed. 


it is interesting to notice that the path in the colour-colour 
plane described by the 100 realizations follows a vector that 
is almost parallel to lines of constant age. At 600 Myr the 
fluctuations cause an almost negligible uncertainty on the 
age, while the effect is more important at 10 Gyr, due to 
the lower sensitivity of {B — K) to age in this age range. 

These fluctuations can clearly be seen in Figs. 4 and 5, 
which show the colours of bright star clusters in the Large 
Magellanic Cloud (LMC) and M31 globular clusters, super¬ 
imposed on our models. The B photometry for the LMC 
clusters is taken from van den Bergh (1981), the J and K 
photometry (plus reddening estimates) is from Persson et 
al. (1983); the age classification comes from van den Bergh 
(1981), and is based on integrated UBV colours. Figure 4 
shows that the young LMC clusters scatter along the tracks 
predicted for the early colour evolution of clusters, and illus¬ 
trates that at these ages (<300 Myr) the colours predicted 
by the models are completely degenerate. The oldest clusters 
lie fairly close to the colour predicted for an age of 14 Gyr 
and [Fe/H] = —1.79, but with a large scatter very similar in 
size and direction to that shown in the simulated old clusters 
in Fig. 3. 

In case of the M31 globular clusters we considered the 
sample recently observed by Rich et al. (2005). The bright 
resolved stars of this cluster sample provide an estimate of 
the individual photometric metallicities and reddenings; the 
integrated colours are taken from Galleti et al. (2004). The 
metallicity range spanned by these M31 clusters is approxi¬ 
mately enclosed by the two metallicities [Fe/H]=—1.79 and 
—0.66 used in our simulations of 10 Gyr old populations. 
Figure 5 shows clearly how the observed colour distribution 



(J-K) 

Figure 5. Observed colours of M31 globular clusters overlaid on 
our model predictions for [Fe/H] between —1.79 and +0.06. Filled 
triangles denote clusters with [Fe/H]< —1.60, while filled circles 
display clusters with higher metallicities. 


follows to a good approximation the distribution of points 
in Fig. 3 for the two old populations. 

As an aside, we have also computed Monte-Carlo sim¬ 
ulations for the {V — K) colour fluctuations, employing the 
same cluster masses and metallicities as an analogous test 
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performed by Bruzual & Chariot (2003), and recovered es¬ 
sentially the same results as these authors. 

To conclude, Figs. 4 and 5 generally show strong encour¬ 
agement for the predictions of our modelling; but it is also 
clear that this statistical fluctuation effect seriously limits 
the use that can be made of stellar clusters in testing these 
predictions. It should be emphasised that these statistical 
fluctuations of the integrated colours are only significant for 
stellar systems up to the sizes of large globular clusters, and 
are completely negligible for galaxies of lO^M© and larger 
(i.e. all galaxies considered in the current paper). 


2.2 Anomalous Horizontal Branch colours 

In general, colours and spectral features involving the blue 
part of the spectrum can be seriously affected by the pres¬ 
ence of blue Horizontal Branch (HB) stars not accounted 
for in the theoretical calibration. Such stars have high tem¬ 
peratures and can add an appreciable contribution to the 
flux e.g. in the B filter and Balmer lines. The models we 
used in our analysis have been computed with a standard 
Reimers (1975) mass loss law along the RGB, with the free 
parameter rj = 0.2. This choice allows one to reproduce the 
mean colour of the HB in Galactic globular clusters of var¬ 
ious metallicities, that qualitatively corresponds to a shift 
from red to blue - for a mean value of the age of ^ 12-13 
Gyr - when the metallicity decreases. However, as has been 
well known for many years (e.g. Sandage & Wildey 1967) 
some Galactic globular clusters of intermediate metallicity 
(e.g. [Fe/H] between ~ —1.6 and ^ —1.1) like the pairs M3- 
M13, NGG362-NGC288, show significantly different mean 
colours (M13 and NGC 288 having a much bluer HB with 
respect to their counterparts, that display an HB colour in 
line with other clusters of the same metallicity) in spite of 
having apparently the same age and the same [Fe/H]. This 
occurrence is the so-called second parameter phenomenon. 
Indeed, one explanation for this different morphology would 
be a higher age for the blue HB clusters (lower evolving 
mass hence bluer HB colours) but this appears not to be 
confirmed by their TO brightness. 

The effect of an unaccounted-for blue HB on the {B — 
K)-{J — K) integrated colours has been evaluated with the 
following numerical test. For ages of 10 and 14 Gyr, and 
[Fe/H]=—1.27 we have computed isochrones using a more 
efficient mass loss along the RGB (the value of rj has been 
doubled) that produced an HB much bluer than the case of 
our standard colours discussed above. The integrated colours 
for these blue HB populations have then been computed ho¬ 
mogeneously with our standard case, and are displayed in 
Fig. 6. As expected, the effect on the inferred age is signif¬ 
icant; these blue HB populations appear to be ^ 6 Gyr old 
when the age is evaluated using our standard calibration, 
instead of 10 and 14 Gyr. The metallicity estimate is only 
slightly affected, because of a small contribution of the HB 
to the J and K colours. This quantitative estimate is typi¬ 
cal for all populations with ages greater than ^10 Gyr, i.e. 
the age range where this effect is found at work in Galactic 
clusters. 



(J-K) 

Figure 6. Model colours, with added points (triangles) show¬ 
ing the effect of, respectively, a bluer Horizontal Branch (for a 
population with t=10 Gyr and 14 Gyr, [Fe/H]=—1.27) an a- 
enhanced metal mixture (for intermediate-old populations with 
[M/H] =—0.66 and 0.06), and a variation of the IMF exponent 
(for two representative populations with [Fe/H]=0.06 and ages of 
600 Myr and 14 Gyr, respectively). Filled circles mark the inte¬ 
grated colours predicted by our standard calibration for the age- 
metallicity combinations used in our experiments. The solid lines 
connect the standard colours with those obtained after account¬ 
ing for a bluer Horizontal Branch, an a element enhancement, or 
a change of the IMF exponent. 


2.3 The effect of an a-enhanced metal distribution 

The scaled solar metal distribution, used in almost all stellar 
evolution computations and adopted in our standard mod¬ 
els, is not universal. Probably the most relevant example of 
this phenomenon is the fact that the initial chemical com¬ 
position of both the old field and globular cluster stars in 
the Galactic halo displays a ratio [a/Fe]~ 0.3 — 0.5 (Carney 
1996), where with a we denote the so-called a elements (O, 
Ne, Mg, Si, S, Ca, Ti). Also elliptical galaxies (e.g. Worthey, 
Faber & Gonzalez 1992) and metal poor stars in other Local 
Group galaxies appear to have been formed with an initial 
[Q;/Fe]>0, although with different degrees of enhancement. 

We have studied the effect of an a-enhanced metal mix¬ 
ture on the integrated {B — K) and (J — K) colours by 
performing the following test. The [a/Fe]=0.4 models by 
Cassis! et al. (2004) extended along the AGB TP phase 
have been employed to determine the integrated colours for 
the following age-metallicity combinations: [M/H] =—0.66, 
age 3 and 10 Gyr; [M/H]=0.06, age 10 Gyr. We re¬ 
call that in the case of a scaled solar metal mixture 
[M/H]=[Fe/H], whereas in case of this adopted a-enhanced 
mixture [Fe/H] = [M/H]—0.35. The results of this test are dis¬ 
played in Fig. 6. The age inferred for the a-enhanced pop¬ 
ulations when using our standard calibration is the correct 
one for the 10 Gyr population, whereas there is an age over¬ 
estimate of about 1-2 Gyr for the 3 Gyr case. The [a/Fe]>0 
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(J — K) colours are shifted to the blue with respect to the 
scaled solar ones, and their values are within ~0.1 dex of the 
(J — K) corresponding to [Fe/H] = [M/H]—0.35 and a scaled 
solar metal mixture. We have found these same results con¬ 
cerning the recovered age and [Fe/H] when we performed 
the same comparisons at both lower and larger metallicities. 
Overall, the combined effect of an a-element enhancement 
on both the theoretical models (i.e., evolutionary timescales 
and the luminosity-Te// plane) and the bolometric correc¬ 
tions to the B, J and K bands, makes this colour combina¬ 
tion a good tracer of the [Fe/H] ratio, even in the presence 
of a ratio [a/Fe]>0. 

To summarize, a scaled solar calibration of the (B — K)- 
(J — K) diagram is able to recover the correct age and ap¬ 
proximately the correct [Fe/H] of an a-enhanced population 
with [a/Fe] typical of the Galactic halo. In case of a younger 
3 Gyr old population, the inferred age tends to be overesti¬ 
mated by 1-2 Gyr. Similar offsets are found over the entire 
range of metallicities investigated in this paper. 

2.4 The effect of the IMF exponent 

To estimate the effect of possible departures from the 
Salpeter IMF, we have computed the integrated colours of 
solar metallicity SSPs using different exponents of the IMF. 
The results for two specific ages (14 Gyr and 600 Myr, re¬ 
spectively) are shown in Fig. 6. As a general rule, a more 
negative value of the exponent (MS dominated population) 
shifts the integrated colours towards bluer values, whereas 
the opposite is true in case of a less negative exponent (giant 
dominated population). 

The reason for this behaviour is easily understood by 
recalling that MS objects less massive than ^ Q.5Mq do 
not contribute appreciably to the total integrated flux in 
the our photometric bands of choice, even if we decrease 
the exponent of the IMF significantly. For example, if we 
set the exponent to —3, excluding the contribution of stars 
with M < O.5M0 decreases the (J — K) colour of a 10 Gyr, 
[Fe/H]=0.06 population by only 0.012 mag and the {B — K) 
colour by 0.08 mag. An MS dominated SSP will have a larger 
contribution of MS stars to the integrated flux, but MS ob¬ 
jects in the relevant mass range (M > 0.5M©) are typically 
bluer than bright RGB and AGB objects (that still pro¬ 
vide the largest contribution to the J and K fluxes) and 
thus make the integrated colours slightly bluer compared 
to the case of the Salpeter exponent. The opposite is true 
in a giant dominated SSP. The magnitude of this effect in¬ 
creases with decreasing age, because the extension of the MS 
reaches brighter magnitudes and bluer colours when the age 
decreases. 


2.5 Composite stellar populations 

The colour-colour diagram of Fig. 1 refers to simple stellar 
populations, made of coeval stars born with the same chem¬ 
ical composition. In the case of a complex star formation 
history, the age and metallicity obtained from this diagram 
will reflect some mean value, that is determined by the num¬ 
ber and luminosity of the objects belonging to various stellar 
generations formed during the system’s lifetime. 

The age sensitive (B—K) colours are obviously very sen¬ 


sitive to the presence of young populations, whose Main Se¬ 
quence is very bright in the B band. The metallicity-sensitive 
{J — K) colour is strongly affected by intermediate age pop¬ 
ulations (300 - 600 Myr old) that show a very bright AGB. 

We found that, keeping the metallicity constant, the age 
inferred from (B — K) is always biased towards the youngest 
stellar component. In the case, for example, of two star for¬ 
mation bursts 40 Myr and 10 Gyr old, the {B — K) is com¬ 
pletely dominated by the youngest population, even if this 
contributes to the total system mass by only ~ 3-5%. This is 
very similar to the conclusions reached by Tantalo & Chiosi 
(2004) in a theoretical study of the age-determinations us¬ 
ing H/l and [MgFe] line indices. Their models showed that 
a burst of star formation involving 2% of the stellar mass of 
an otherwise old population would significantly perturb the 
derived age to young values for a period of order 2 Gyr after 
the onset of the burst. 

A powerful test for the adequacy and usefulness of the 
theoretical {B — K) — {J — K) calibration discussed in this 
paper, is provided by the comparison of our model predic¬ 
tions with integrated colours of Local Group galaxies, whose 
star formation histories (star formation rate - SFR - plus 
age-metallicity relation - AMR) are independently estimated 
from their resolved stellar populations. 

2MASS J and K total integrated magnitudes for the 
Local Group galaxies NGG 6822 (dwarf Irregular) and 
NGG 185 (dwarf Elliptical) have been taken from Jarrett 
et al. (2003), while the integrated B magnitudes have been 
taken from the Third Reference Gatalogue of Bright Galax¬ 
ies (de Vaucouleurs et al. 1991). Estimates of the star forma¬ 
tion histories of these two galaxies, using Golour-Magnitude- 
Diagrams of their resolved stellar populations, are available 
in the literature (i.e., Gallart et al. 1996, Wyder 2001 for 
NGG 6822; Martinez-Delgado, Aparicio & Gallart 1999 for 
NGG 185). Observed (J — K) colours, and K magnitudes 
have been transformed to the Bessell & Brett (1988) system 
using the transformations by Garpenter (2001). Extinction 
corrections from Schlegel, Finkbeiner & Davis (1998) have 
been applied, and no correction for possible internal extinc¬ 
tion has been considered. 

Figure 7 displays the integrated colours of these two 
galaxies compared to our model grid. NGG 6822 appears 
to be young, located in the region of full degeneracy be¬ 
tween age and metallicity effects. In light of our previous 
discussion about composite stellar populations, the position 
of NGG 6822 on our colour-colour diagram is fully consis¬ 
tent with its estimated star formation history, that displays 
a sizable star formation rate at young ages. NGG 185 ap¬ 
pears to be older, with a mean age slightly above 3 Gyr and 
a mean metallicity of the order of [Fe/H]~ —1.8. This is 
again consistent with the star formation rate inferred from 
its resolved stellar population, dominated by ages above ~ 
8 Gyr, and a much smaller star formation activity at young 
ages. 

A more quantitative test has been performed in case of 
NGG 6822 (an analogous quantitative estimate for NGG 185 
could not be performed, because of the lack of a firm esti¬ 
mate of the AMR for this galaxy). Wyder (2001) provides 
the SFR and AMR for three fields within the central bar, 
and two outer fields. The star formation history of the cen¬ 
tral bar fields - located in a region that contributes most to 
the total integrated magnitude of the galaxy - has been em- 
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Figure 7. Observed colours of the Local Group galaxies NGC 185 
and NGC 6822 overlaid on our model predictions. Crossed circles 
denote theoretical integrated colours obtained from the SFR and 
AMR of NGC 6822, as estimated by Wyder (2001 — see text for 
details). 

ployed to determine their expected integrated colours. These 
theoretical predictions, displayed in Fig. 7, show a very good 
agreement with the observed integrated colour of the whole 
galaxy. 


3 APPLICATION OF THIS METHOD TO 

EARLY-TYPE GALAXIES 

In this section we will use the model colours to determine 
properties of the early-type galaxy population in the nearby 
Virgo cluster and the general field. The 2MASS database 
contains JHK photometry for most of the giant galaxies in 
Virgo, and a smaller but useful number of the brighter dwarf 
galaxies, so the main problem is finding matching R-band 
photometry. Our initial analysis made use of Bt extrapo¬ 
lated total magnitudes from the Virgo Cluster Catalogue 
(Binggeli, Sandage & Tammann 1985), which we combined 
with the 2MASS K^xt magnitudes to derive total B — K 
magnitudes (see e.g. Mannucci et al. 2005). However, we 
found systematic offsets between the colours thus obtained 
and those calculated from aperture photometry of the same 
galaxies (Poulain 1988, Poulain& Nieto 1994), in the sense 
that the extrapolated magnitudes gave bluer colours than 
aperture measurements, typically by 0.2-0.3 mag., so we 
restricted the analysis to galaxies with matched aperture 
photometry in all three passbands {B, J and K). 

One extremely useful database is that of Michard (2005) 
who studies the optical/near-infrared colours of 94 elliptical 
galaxies out to a distance modulus (m — M)o=33.52 and 
with total integrated B magnitude Bt <18.8, using 2MASS 
near-IR images and optical aperture photometry from the 
literature. He quotes colours within an effective aperture Ae; 



Figure 8. Integrated colours of a sample of nearby elliptical 
galaxies overlaid onto the theoretical calibration. Selected ages 
and metallicities are marked; the typical observational (Icr) er¬ 
rors and the reddening vector are also shown. Filled circles show 
the galaxies belonging to the Virgo cluster. 


in every case this is interpolated from the published multi¬ 
aperture photometry in exactly the same way for the optical 
and near-infrared photometry, so no systematic colour dif¬ 
ferences should result from this procedure. Note that the B 
photometry we use, which is originally from Poulain(1988) 
and Poulain & Nieto (1994), is not mentioned in the main 
text of Michard’s paper, but is available in the associated on¬ 
line data. The integrated magnitudes have been dereddened 
and fc-corrected as described in Michard (2005), and the in¬ 
frared magnitudes and colours transformed to the standard 
Bessell & Brett (1988) system following Carpenter (2001). 

Figure 8 displays the colours of the 94 galaxies overlaid 
onto the theoretical calibration. The filled symbols show 16 
objects belonging to the Virgo cluster, the highest density 
environment in the region of space sampled by these data. 
The distribution of about 99% of the galaxies in the sample 
follows nicely a mean linear relationship between {B — K) 
and (J — K), corresponding to an average age of 8 Gyr 
and a range in [Fe/H] values between -1-0.4 and about the 
solar value. Superimposed on this mean relationship is an 
age spread that spans the range between ^3 and ~14 Gyr. 
As an example, in the Virgo cluster NGG 4478 has a lu¬ 
minosity weighted age of about 3 Gyr, whereas NGC 4374 
is about 10 Gyr old. These age and metallicity spreads ap¬ 
pear to be real, and are not due to the observational errors, 
which are quite small in comparison. We note that the use 
of the models of Girardi et al. (2002) would result in some¬ 
what higher mean metallicities and younger mean ages in 
the colour range covered by this galaxy sample. 

The distributions of points corresponding to the Virgo 
galaxies and the rest of the objects are statistically the same, 
as verified with a 2-dimensional Kolmogorov-Smirnov (KS) 
test (Press et al. 1992). Thus we find no evidence for envi- 
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Figure 9. Radial colour gradients for selected galaxies. The upper 
plot shows 8 bright Virgo cluster elliptical galaxies, which are 
identified in the text. In each case, the range of aperture colours is 
indicated by the two points plotted for each galaxy, with the right- 
hand point indicating the colours in the smallest aperture and 
the left hand point those for the largest aperture. The lower plot 
shows similar results for the 5 dwarf galaxies, also identified in the 
text. Again, the two points indicate the range of colours found, 
with the reddest colours being generally found for the smallest 
apertures. The clear exception to this rule is NGC 4640, for which 
all available aperture colours are plotted explicitly. In this case, 
the smallest aperture (10” diameter) has the reddest B—K colour, 
and this colour gets monotonically bluer to the largest aperture 
(50”). However, there is no overall trend in J — A for NGC 4640. 


ronmental effects on the luminosity-weighted mean ages and 
metallicities of bright ellipticals in the sample of Michard 
(2005). 

In general, the integrated {B — K) — {J — K) colours 
of this sample show that, if ellipticals formed in a single 
star formation episode, their epoch of formation spans a 
large range of redshifts. If multiple generations are present 
together with an old component formed soon after the Big 
Bang, the age range derived in this analysis can be explained 
in terms of different fractions of intermediate age/young 
populations, which bias the luminosity-weighted ages to¬ 
wards lower values. The presence of very young populations 
(ages below 100 Myr) at a level larger than ^1-2% (by mass) 
can be excluded by the results of this colour-colour analysis. 


4 AGE AND METALLICITY GRADIENTS 

WITHIN GALAXIES 

One complication in comparing results from integrated 
colours with literature studies based on line indices is that 
the latter are generally based on spectroscopic measure¬ 
ments which are only available for central regions of the 
galaxies concerned. Thus it is important to check for sys¬ 
tematic gradients in derived properties of stellar populations 
which may affect such comparisons. One advantage of the 
methods presented here is that it is simple to get spatially- 
resolved information from array imaging, and hence to com¬ 
pare central and global stellar populations. 

A subsample of 8 Virgo galaxies with aperture pho¬ 
tometry from Poulain (1988), Poulain & Nieto (1994) and 
2MASS is displayed in Fig. 9, upper frame; the sample con¬ 
tains NGC 4261 (open circles), NGC 4365 (filled circles), 
NGC 4374 (open pentagons), NGC 4472 (open squares), 
NGC 4473 (filled squares), NGC 4486 (filled pentagons), 
NGC 4621 (filled triangles) and NGC 4636 (open triangles). 
The aperture sizes (diameters) used by these authors were 
21"9, 31”2, 43”4, 60”9 and 86^6 (although in the case of 
NGC 4473 the largest aperture is 60”9). However, for every 
one of these galaxies we found a monotonic trend in colours 
with aperture, with the reddest colours being found for the 
smallest aperture, so we only plot symbols for the smallest 
and largest apertures. The typical error bars associated to 
the colours through an individual aperture are displayed in 
the upper left corner. 

The general trend is of redder colours for decreasing 
aperture, a well known result (see, e.g. Peletier et al. 1990). 
The two galaxies with the reddest (J — K) colours are 
NGC 4261 and NGC 4486, both of which are radio galaxies 
with strong nuclear emission, and so the colours are not re¬ 
liable indicators of the underlying stellar population. Inter¬ 
estingly, for NGC 4374, NGC 4473, NGC 4636 and, to some 
extent, also NGC 4621 (if one takes into account the size 
of the error bars on the colours of each individual aperture) 
the colour gradients correspond to a sequence of increasing 
metallicity at approximately constant age. The slope of these 
radial colour gradients in the {B—K) — {J—K) plane is there¬ 
fore similar to the mean slope of the total integrated colours 
for the whole sample of galaxies, a result in agreement with 
Peletier et al. (1990) who employed the {U — V) — {V — K) 
plane. Approximately constant luminosity-weighted ages at 
different radii are also found in relatively young looking 
galaxies like NGC 4636 and NGC 4473 (luminosity weighted 
ages of ~6 Gyr and ~3 Gyr respectively); this implies that 
any intermediate age/young population that eventually co¬ 
exists with an old («10-14 Gyr) first stellar generation must 
be evenly spread throughout the parent galaxy. 

The situation for NGC 4365, and to some extent 
NGC 4472, is different. Moving from the largest aperture 
towards the smallest one, the {B — K) colour stays approx¬ 
imately constant, whereas (J — K) increases. This implies 
that the luminosity-weighted age decreases from ~ 10 Gyr 
to below 6 Gyr, and [Fe/H] increases from ~ 0.06 to above 
0.40. Internal extinction should not be responsible for these 
different colour gradients (nor for the colour gradients of all 
other objects displayed in the figure) given that these are all 
early-type galaxies. Only a progressively bluer morphology 
(which would have to be bluer than the underlying theo- 
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retical calibration) of the HB with decreasing radius would 
be able to keep the NGC 4365 {B — K) colour constant for 
increasing (J — K) and metallicity. 

The lower frame in Fig. 9 shows similar results for 5 
dwarf galaxies in the Virgo cluster. These are the only Virgo 
cluster galaxies classified as dwarfs by Binggeli et al. (1985) 
which have 2MASS photometry with errors in J — K < 
0.08 mag., and for which B-band photometry is available 
from the survey of Sabatini et al. (2003). The latter study 
covered two strips of total area 25 deg^ within the Virgo clus¬ 
ter. The galaxies plotted are UGC 7436, type dE5, plotted 
as filled squares; IC 781, type dS0,N, plotted as Hlled trian¬ 
gles; IG 3292, type dSO, plotted as filled circles; UGG 7399A, 
type dE5,N, plotted as open circles; and NGG 4640, type 
dS0,N, plotted as open squares. Again, only the smallest 
(10") and largest (40-60") apertures are plotted for most 
of these; where this is done, the smallest apertures always 
have the reddest colours. Only in the case of NGC 4640 have 
all the aperture colours been plotted separately. For some 
reason this galaxy exhibits very different colour gradients 
from the other dwarfs. In this case the B — K colours red¬ 
den monotonically from the smallest aperture (lO") to the 
largest aperture (50"), while the J — K colour does a pecu¬ 
liar loop to redder colours and back again. Apart from this 
galaxy, the overall colour gradients appear to be dominated 
by increasing metallicity towards the centres of these galax¬ 
ies, with any age gradients being small, and those which are 
seen going in both directions. However, the most interesting 
result from this small number of dwarf galaxies is not in the 
colour gradients, but in the overall position of these dwarfs 
compared to the more luminous galaxies seen both in the up¬ 
per frame of Fig. 9 and in Fig. 8. The dwarfs are shifted with 
respect to the brighter galaxies, both towards lower metal- 
licities and, more interestingly, towards lower luminosity- 
weighted mean ages, with a range for the dwarfs of 1-6 Gyr. 
It is hard to draw any general conclusions from such a small 
sample, but the difference is quite striking, particularly given 
that we are studying the very brightest end of the dwarf 
galaxy luminosity function, because of the difficulty in de¬ 
tecting lower luminosity and hence lower-surface-brightness 
galaxies with the relatively shallow 2MASS data. 

Thus we do see significant gradients in colours as a func¬ 
tion of aperture size, which needs to be borne in mind when 
comparing integrated properties using colours with nuclear 
properties using line indices, particularly since spectroscopic 
apertures will tend to be much smaller even than the small¬ 
est apertures used here. 


5 DISCUSSION AND CONCLUSIONS 

We will summarise briefly our conclusions on the nature 
of the stellar populations of the various galaxy samples as 
revealed by their optical/near-infrared colours. It is useful 
to remember at this point that the ages and metallicities 
derived are luminosity-weighted mean values, and we have 
shown that, as for spectroscopic indices, the results can be 
heavily skewed by the presence of even a small (a few per 
cent by mass) young stellar population in an otherwise old 
galaxy. 

For the bright elliptical galaxy sample of Michard 
(2005), we find that the integrated colours within the ef¬ 


fective radius for these galaxies reveal stellar populations 
with mean ages averaging about 8 Gyr, but with a real scat¬ 
ter, larger than the measurement errors, covering the range 
from 3-14 Gyr. Mean [Fe/H] values are just above solar, 
but again there is a real scatter, from just below solar to 
just above [Fe/H]=0.4. The Virgo cluster galaxies contained 
within Michard’s sample show no statistically significant dif¬ 
ferences in their integrated parameters than the sample as 
a whole, so there is no evidence for strong environmental 
effects on either mean stellar ages or metallicities. 

A more detailed study of population gradients within 8 
Virgo cluster giant galaxies and 5 Virgo cluster dwarf galax¬ 
ies reveals a strong trend for central regions of all galaxies 
to have higher metallicities than the galaxies as a whole. 
Evidence for age gradients is much more marginal, but one 
giant galaxy shows evidence for a younger central popula¬ 
tion, while the dwarfs show weak evidence for age gradients 
in both directions. 

Gonsidering the 5 Virgo dwarf galaxies together, we hnd 
that they have younger mean stellar ages than their more lu¬ 
minous counterparts. Glearly, this is too small a sample to be 
drawing general results about the dwarf galaxy population 
as a whole, but it should be noted that such a result would be 
consistent with recent observational results by Caldwell et 
al. (2003), using spectral indices of low-luminosity galaxies 
in Virgo and lower density environments, and Kauffmann 
et al. (2003), who perform a statistical study of the star 
formation histories of galaxies from the Sloan Digital Sky 
Survey. 

There are two conclusions we can make from this pre¬ 
liminary study. Firstly, if we assume that all these galax¬ 
ies contain a simple stellar population (SSP), then quite 
contrary to the hierarchical model of galaxy formation the 
smallest galaxies have not formed first at redshifts of 6-10 
(ages greater than 12 Gyr for standard ACDM) but, this 
does ht in with the ‘downsizing’ scenario discussed in the 
introduction. Secondly, if these galaxies do not contain a 
SSP but do have a small fraction of young stars then the 
feedback mechanisms discussed in the introduction are in¬ 
efficient - these galaxies have retained at least some gas to 
enable them to continue to form stars over long periods of 
time (see also the model by Kobayashi, 2005, which has this 
property). This second interpretation fits in very well with 
observations of the resolved stellar populations found in Lo¬ 
cal Group dwarf galaxies (Grebel & Stetson 1999; Smecker- 
Hane et al. 1996; Gallart et al. 1999; Mighell & Rich 1996) 
and discussed in the introduction. 

Unfortunately, it is only possible to study the brightest 
of the Virgo dwarf elliptical galaxies using 2MASS photom¬ 
etry. However, with new data becoming available, e.g. from 
the UKIDSS project, it should be possible to derive the re¬ 
quired colours for some hundreds of Virgo cluster dwarfs, 
pushing significantly further down the luminosity function. 
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